Relation between phase space coverage and entanglement for spin-1/2 systems 
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For systems of two and three spins 1/2 it is known that the second moment of the Husimi 
function can be related to entanglement properties of the corresponding states. Here, we generalize 
this relation to an arbitrary number of spins in a pure state. It is shown that the second moment 
of the Husimi function can be expressed in terms of the lengths of the concurrence vectors for all 
possible partitions of the A'^-spin system in two subsystems. This relation implies that the phase 
space distribution of an entangled state is less localized than that of a non-entangled state. As 
an example, the second moment of the Husimi function is analyzed for an Ising chain subject to a 
magnetic field perpendicular to the chain axis. 
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I. INTRODUCTION 

In the last few years there has been growing activ- 
ity in the study of the behavior of spin chains viewed 
from a quantum information perspective. The relations 
between condensed matter physics and quantum infor- 
mation are twofold in this case. On the one hand, spin 
chains can provide a tool for quantum communication 
[l|, On the other hand, the concept of entanglement 
has been employed to study spin systems, in particular at 
a quantum phase transition [3,|3jS[6'|. Furthermore, the 
concept of matrix product states has led to new insights 
into the density matrix renormalization group algorithm 
(DMRG) with which the ground state properties of spin 
systems can be determined [7, 8]. 

Recently, A. Sugita @ has pointed out for two spins 
1/2 a relation between an entanglement measure, the so- 
called concurrence [l(]| , and a property of the phase-space 
representation of the state of the two spins. More specif- 
ically, this relation involves the second moment of the 
Husimi function, a positive definite phase space distribu- 
tion, [ill , [l^ This quantity can be viewed as an inverse 
participation ratio in phase space and measures to which 
extent the phase space is covered by the Husimi function. 
It turns out that the more the state extends over phase 
space, the more the state is entangled. In particular, the 
Husimi function of a factorizable state minimally covers 
the phase space. A first impression of this difference be- 
tween the phase space representations of entangled and 
non-entangled states can be obtained from Fig.[T] In this 
figure, the basic structure of the Husimi function is vi- 
sualized by full lines where the maxima of the Husimi 
function are located. At dashed lines and on gray planes 
the Husimi function vanishes. The comparison of a fac- 
torizing state on the left and a maximally entangled state 
on the right indicates that in the latter case the exten- 
sion of the Husimi function is larger. The figure will be 
explained in more detail in Sec.lTTlbelow where these first 
observations will be made more precise. 

For three spins an expression for the second moment 
of the Husimi function was given by Sugita in terms of 
the concurrence between two of the three spins and the 
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FIG. 1: Basic structure of the Husimi function of a factorizing 
state (left) and a Bell state (right). The full lines indicate 
the positions of the maxima of the Husimi function. Along 
the dashed lines and on the gray planes the Husimi function 
vanishes. The entangled state has a more complex phase space 
structure than the non-entangled state. For details see Sec. [Hi 



so-called 3-tangle As the Husimi function can be de- 
fined for an arbitrary number of spins, it is interesting to 
determine whether and if yes how its second moment is 
related to entanglement measures for an arbitrary num- 
ber of spins. 

Another motivation for this study arises from the fact 
that recently phase space methods have been employed 
to analyze condensed matter systems. In particular the 
occurrence of a metal-insulator transition in disordered 
and quasiperiodic systems has been analyzed by means 
of the inverse participation ratio in phase space [3, [IB] ■ 
It appears natural to extend these studies to interact- 
ing systems. Unfortunately, already the numerical treat- 
ment of two particles moving in one dimension is quite 
demanding. Here, spin systems like the one which we will 
discuss in Sec. IIVI present significant advantages. For the 
interpretation of the phase space properties, it is again 
interesting to know the relation between the second mo- 
ment of the Husimi function and the entanglement of the 
state under consideration. 

The phase space approach has occasionally been used 
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in the context of quantum information in the past. The 
state and time evolution of a quantum computer have 
been described by means of a discrete Wigner function 
[TgI [iTj . More closely related to the present work is the 
investigation of the entanglement in bipartite systems on 
the basis of the Wehrl entropy where each subsystem is 
described in terms of a sufficiently large spin [l8|. Fur- 
thermore, the relation between nonlinear bipartite sys- 
tems and the classical phase space dynamics has been 
considered In contrast, here we consider the phase 
space properties of an arbitrary number of distinguish- 
able spins by means of a corresponding number of spin- 
1/2 coherent states. A single number extracted from such 
a phase space representation can at best describe entan- 
glement in a very global sense. Therefore, our intention 
is not to introduce a new quantity to measure entangle- 
ment but rather to provide a link between phase space 
properties and entanglement. 

We start in Sec. |TT] by reviewing the phase space rep- 
resentation for spin- 1/2 systems and presenting Sugita's 
results in a form suitable for the ensuing discussion. In 
Sec. mil the second moment of the Husimi function will 
be expressed in terms of projectors acting on the Hilbert 
space and an auxiliary copy. In this way, we make con- 
nection to the work by Mintert et al. [2^ which allows us 
to obtain a relation to concurrences. Finally, in Sec. IIVI 
we discuss as an illustrative example the Ising spin chain 
in presence of a magnetic field perpendicular to the chain 
axis. 



II. PHASE SPACE FOR SPINS 

We will restrict our discussion to spin- 1/2 systems 
where the spin-coherent states, the analogue of the co- 
herent states for the harmonic oscillator, are given by 

MM 

|^,V:,)=cos(^0|O)+sin(^0e'^|l) (1) 

so that each point on the Bloch sphere characterized by 
the two angles and ip represents a spin-coherent state. 
We employ the notation common in quantum informa- 
tion, where |0) and |1) correspond to the eigenstates of 
the Pauli matrix az with eigenvalue -1-1 and —1, respec- 
tively. A coherent state for a system consisting of N 
distinguishable spins can be expressed as product of co- 
herent states for each of the spins 

N 

(2) 

■i=l 

The spin-coherent states enable us to define a pos- 
itively definite phase space distribution 

H{fi) = {fi\p\fi) (3) 

which is the spin analogue of the Husimi- or Q-function 
familiar from the harmonic oscillator. [ll|, [l2| In ([31), 



p denotes the density matrix of the state for which the 
phase space distribution is determined. 

In order to quantify the extension of a state in phase 
space, we introduce the second moment of the Husimi 
function 

V = 3^ J dpHifif (4) 

with the Haar measure dfi = J^^j^ sin(?9i)d?9i(i(pi/47r. 
The prefactor is chosen such that a separable state leads 
to = 1. Up to a factor (3/2)^, V is the inverse par- 
ticipation ratio in phase space. Its inverse measures the 
extension of the Husimi function in phase space. We 
note that the second moment of the Husimi function cor- 
responds to the first nontrivial term in the expansion of 
the Wehrl entropy J dpniii)\n[niii)]. ^ 

To get a feeling for the physical content of the Husimi 
function Ti and its second moment V, we review results 
obtained for systems containing two or three distinguish- 
able spins. For a pure two-spin state 

IV-) =a|00)+6|01)+c|10)+d|ll) (5) 

with |ap + |6p + |cp + \d\'^ = 1 the second moment of the 
Husimi function is obtained as 

r = l-\ad-bc\^ . (6) 

For separable states, one has ad — be and therefore V = 
1. On the other hand, the minimal value of the second 
moment of the Husimi function for systems consisting of 
two spins is obtained for Bell states with V — 3/4. These 
results indicate the existence of a relation between this 
phase space quantity and the amount of entanglement. 

In order to illustrate the difference between the two 
cases, it is useful to first consider the Husimi function 
which, for two spins 1/2, in general is a function of the 
four angles 'di,'&2,Vi, and (p2- For the factorizing state 
lip) — 1 00), one obtains the Husimi function 

H= i(H-cos(i9i))(l + cos(z92)) (7) 
while for the BeU state I?/;) = (|00) + |ll))/\/2 one finds 

H = i(l+cos(t?i) cos(t?2)+sin(t?i) sin(i92) cos{ipi+ip2)) ■ 

(8) 

The fact that these Husimi functions depend only on 
three independent variables allows us to represent their 
structure in Fig. [T] For the factorizing state |00), the 
maximum depicted by a full line lies at i9i = = 
as should be expected on the basis of ([1]). The gray ar- 
eas delimited by the dashed lines indicate the planes at 
-di — TT and i?2 = tt where the Husimi function vanishes. 
Similarly, for the Bell state in Fig. [Jd, the maxima at 
"i^i = 1^2 = and TT indicate the presence of the states |00) 
and 1 11). The relative phase between the two states is en- 
coded in the position of the bridges along -ffi = ??2- The 



3 



dashed lines again indicate zeroes of the Husimi function. 
We remark that the independence of the Husimi function 
on (pi — ip2 is specific to superpositions of the states |00) 
and 111). 

The second moment of the Husimi function ([6|) for two 
spins is related to the concurrence [l^l 



C = 



(9) 



where the star denotes the complex conjugate in the 
eigenbasis of CTz. An alternative expression, which will 
be useful in the following, can be obtained by introduc- 
ing an auxiliary Hilbert space with a copy of the state 
We denote the state in both Hilbert spaces as a 
column vector |^). Then the concurrence becomes 



C = 
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(10) 



where the first bracket operates in the Hilbert space of 
the first spin while the second bracket operates in the 
Hilbert space of the second spin. A more general discus- 
sion of expressing the concurrence in terms of projectors 
can be found in Ref. [l^- Beyond these formal consider- 
ations, such an auxiliary Hilbert space has very recently 
been employed to directly measure the concurrence of 
two photons. [25j 

For the general two-spin state ([5]), the concurrence is 
given by 



C = 2\b*c* - a*d* 



(11) 



so that by comparison of ([6]) and pT|) one immediately 
obtains the relation 



V = 1 
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(12) 



For three spins, the second moment of the Husimi func- 
tion can be expressed in terms of the concurrences be- 
tween two of the three spins A, B, and C and the 3-tangle 
TABC [il as 



^ = 1-J(C1b 



(13) 



For the generalization to an arbitrary number of spins it 
is more suggestive to express this result in terms of the 
concurrences between one spin and the two others as 



^ - 1 ~ g('^l(BC) 



'^B(AC) 



' '^C(AB) 



)• (14) 



We remark here that the cases of two and three spins 
are the simplest in the sense that each partition into two 
subsystems will yield a single concurrence. This is in 
general no longer true for more than three spins, the case 
which we are going to address now. 



III. SECOND MOMENT OF THE HUSIMI 
FUNCTION AS A PROJECTION 

For the following considerations it is convenient to ex- 
press the second moment of the Husimi function ^ in 
terms of projectors onto symmetric and antisymmetric 
subspaces. We start by considering a system consisting 
of a single spin 1/2. The key idea is to express the square 
in the integrand of ^ in terms of a tensor product of the 
spin Hilbert space and an auxiliary copy of this Hilbert 
space. The similarity with the auxiliary Hilbert space 
introduced in ([TO)) already hints at the possibility of a 
general relation between the second moment V and the 
concurrence C. 

In order to distinguish between tensor products of dif- 
ferent spins which we note horizontally, the tensor prod- 
uct between one spin Hilbert space and its auxiliary copy 
will be denoted vertically. The density matrix g = p® p 
refers to the tensor product of the density matrices in 
these two spaces. The second moment of the Husimi 
function can then be written as 



■P = Tr 3 dpQ 



(15) 



By construction, the states ||^) do not contain contribu- 
tions antisymmetric under exchange of the two Hilbert 
spaces. 

Expressing the projectors in terms of the coherent 
states ([1]), one can carry out the integrals over the angles 
d and Lp. It turns out that the second moment of the 
Husimi function can be expressed as 



(16) 



Here, = P, 



Pi' 



(1) 



P, 






(2) 



Pi^^ with 



Pi^ 



p(3) ^ i 

2 



(17) 



projects onto the symmetric cigcnstates of two qubits 
while 



Pn - 1 - P, 



(18) 



projects onto the antisymmetric eigenstate. Expectation 
values like in ([TO]) are always to be understood in the 
extended space containing the original Hilbert space as 
well as a copy. 

It is straightforward to generalize this reasoning to 
more than one qubit because a coherent state according 
to ([2|) is defined as a product of coherent states for each 
qubit. For N qubits, the expression (jl6p then becomes 



(P 



(19) 
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Making use of the decomposition 1 = (P^ + Pa)®^ we 
can write this expression in the more compUcated but 
useful form 



where 



+---+{pr}). 



PH (20) 



The curly braces imply a sum over all different orderings 
of projection operators. Noting that 



P. - P, 



one can demonstrate the relation 

Tr {g{P, - Pa)^^) = Tr(p2) 

which allows us to rewrite ([20]) as 
1 



(21) 



(22) 



P = -(1 + Tr(p2)) 



^pn +{pf''~'^pn +■■■)■ 

(23) 

From it follows that for mixtures all terms in ([^0]) 
will contribute while for pure states the terms with an 
odd number of projectors Pa onto antisymmetric states 
are irrelevant. All these contributions clearly vanish be- 
cause a state |^) is symmetric under exchange of the 
Hilbert space and its auxiliary copy. 

For the further discussion, we will restrict ourselves 
to pure states where the second moment of the Husimi 
function now reads 

V = l- ({Pf <E> PD + {Pf <E> Pn +■••)• 

(24) 

In particular, for at most three qubits only one term in 
the expectation value will contribute, yielding the simple 
relations ([H]) and P^ . 

The expression (j24p depends on a linear combination 
of projectors with equal weight. It represents a special 
case of a class of operators which can be employed to 
define a concurrence. Following Ref. [13, we introduce 
the iV-partite concurrence of a pure state |^) describing 
a system consisting of N spins 



cnW = 21-^/^(2^^ - 2)(m' ~T.^^p1 (25) 

and find for the second moment of the Husimi function 

cnW 



1 



(26) 



In (|25p . Pi is the reduced density matrix of a subsystem 
and the sum runs over the 2^ — 2 subsystems containing 
at most iV — 1 spins. Alternatively, one can make use of 
the relation [lO, [2611 



(27) 



1/2 



c = 



[partitions a 



(28) 



describes the total length of the concurrence vectors for 
all partitions of the system into two subsystems. The 
index a denotes the components of the concurrence vec- 
tor for a given partition. We thus arrive at the relation 
between the second moment of the Husimi function and 
the total length of the concurrence vectors 



V = 1 



(29) 



This relation generalizes the results and (fT^l) which 
are immediately recovered by noting that for two and 
three spins each concurrence vector contains only one 
component and that there exist one and three partitions, 
respectively. 

From (p9|) one can conclude, that entanglement leads 
to a decrease of the second moment of the Husimi func- 
tion and thus to a larger spread of the phase space distri- 
bution. The relevant measure of the entanglement here 
is the total length of the concurrence vectors. 

It is instructive to determine the second moment V for 
two different entangled states, the iV-qubit GHZ and W 
states. For the A^-qubit GHZ state 



I GHZ 



V2 



(loo-.-o) 



one finds 



P(|GHZ^)) 



111 



1 

2^ 



1)) 



(30) 



(31) 



This result makes sense even for only one or two qubits. 
In the first case, one finds P = 1 indicating a "separa- 
ble" state while in the second case the second moment of 
the Husimi function of a Bell state is recovered. As the 
number of qubits increases, V decreases and approaches 
the value of 1/2 in the limit N oo. 
For the A^-qubit W state 

|Wjv) = -^(|100---0) + |010---0) + ---) 



= — [|1) ® |00 • • • 0) + |0) ® vTT^IWn^i)] 
V A 

(32) 

one can derive the recursion relation 

^(IWa^)) = ^^^]^^(|Wiv-i) + ^ . (33) 

With the initial condition P(|Wi)) — 1 one arrives at the 
solution 



P(|Ww)) 



1 
2N 



(34) 
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As for the GHZ state, the case of two qubits reproduces 
the second moment of the Husimi function of a BeU state. 
With an increasing number of qubits, V decreases and 
approaches 1/2 in the hmit N — oo. However, for aU 
N > 2, the second moment of the Husimi function of a 
W state is larger than that of a GHZ state. This reflects 
the fact that the concurrence cpf for a GHZ state is larger 
than that of a W state. [21] Although V does not aUow to 
distinguish the different kinds of entanglement present in 
the two classes of states, a difference is nevertheless visi- 
ble in the phase space structure which is more extended 
for a GHZ state. 

From the results ((3T|) and (p4)) one might infer that 
V possesses a lower bound of 1/2. This is however not 
the case. According to (fT9|) . the second moment of the 
Husimi function for a system consisting of subsystems 
not entangled among each other is given by the product 
of the respective V's of the subsystems. For N pairs of 
spins in a Bell state, one finds r = (3/4)^ which clearly 
goes to zero for iV ^ oo. 

IV. ISING MODEL IN A MAGNETIC FIELD 



fx. 
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FIG. 2: The second moment V of the Husimi function of 
a one-dimensional ferromagnetic Ising model with 8 spins 
and periodic boundary conditions is shown as a function 
of the coupling constant g for various angles between the 
z-axis and the magnetic field taking the values Q/tt = 
0, 0.42, 0.46, 0.48, 0.49, 0.495, 0.4975, 0.4995, and 0.5 from the 
upper to the lower curve. 



In the last few years the connection between entangle- 
ment and quantum phase transitions has been studied 
extensively, for example see Refs. H, |3, S H HI, In 
most cases the concurrence has been used as a measure 
for bipartite entanglement. Taking a different perspec- 
tive we concentrate on the phase space properties of such 
a transition. In the previous section we have seen how 
the concurrence is related to the second moment of the 
Husimi function. It is therefore interesting to study this 
phase space quantity for a model exhibiting a quantum 
phase transition. 

As an example we consider a one-dimensional chain of 
spins 1/2 described by the Hamiltonian 

^ = --^ E + 9 « cos e + < sin e)] . (35) 

i 

The first term for J > gives rise to a ferromagnetic cou- 
pling between neighboring spins while the second term 
arises due to a magnetic field which is oriented at an 
angle O with respect to the 2;-axis. The parameter g de- 
scribes the ratio between the magnetic field strength and 
the interaction strength between two neighboring spins. 
For the case of a transverse magnetic field, Q = 7r/2, this 
model undergoes a quantum phase transition a.t g = 1. 

P 

For 8 = 0, the ground state of ([55)) will be given by 
a factorizing state either with all spins in state |0) or |1) 
depending on the sign of g. Then, V = I independent 
of g. For angles < 8 < 7r/2 and <? <C 1, the spins will 
mostly be in the factorizing state |0 . . . 0). However, as g 
increases, entanglement is build up and the extension of 
the state in phase spaces increases. Correspondingly, V 
will decrease. On the other hand, for g 3> 1 the ferromag- 
netic coupling becomes irrelevant. Then, all spins point 



in the direction of the magnetic field and V should reach 
an asymptotic value of 1. For a transverse magnetic field, 
8 — 7r/2, the ground state in the thermodynamic limit, 
N oo, ior g < 1 will be a GHZ state. According to 
([3T|) . we expect V = 1/2. On the other hand, for g ^ 1, 
the asymptotic value V = I should again be reached. 

In Fig. [21 we present numerical results for a system 
of 8 spins and angles 8 varying from to tt/2. From 
the numerically obtained ground state of ()35p . the sec- 
ond moment of the Husimi function has been determined 
by evaluation of ([T9[) . The lowest curve represents the 
case of the transverse Ising model. For g 1, one finds 
a value for V very close to 1/2 as expected from ()3ip . 
For a finite number of spins, V displays a slight increase 
even below g = I where in the thermodynamic limit V 
is expected to remain at a value of 1/2 before increasing 
for 5 > 1 and reaching V = I asymptotically. Even for 
angles < 8 < 7r/2, the curves in Fig. [2] clearly show 
that entanglement is built up around g = 1. The precur- 
sor of the quantum phase transition thus manifests itself 
also in the phase space properties for angles close to 7r/2. 

We remark that while here we have focussed on the 
ground state properties, the second moment of the 
Husimi function can also be determined at finite tem- 
peratures. Extensions to the antiferromagnetic coupling 
and the two-dimensional Ising model are possible within 
the limits imposed by finite computer resources. (3l| 



V. CONCLUSIONS 

The relation between the second moment of the Husimi 
function and the concurrence derived in for two and 
three spins has been generalized to an arbitrary number 
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of spins. The extension of a state in phase space is thus 
related to its global entanglement properties. Generally, 
entanglement will imply a delocalization of the Husimi 
function of a state. Furthermore, the result (pS)) provides 
a phase space interpretation of the A^-partite concurrence 

The relation between phase space properties and en- 
tanglement has been illustrated by calculating the second 



moment of the Husimi function for the one-dimensional 
Ising model with magnetic field. The precursor of a quan- 
tum phase transition is clearly seen in the presence of a 
transverse magnetic field even for a relatively small num- 
ber of spins. For fields deviating from the transverse di- 
rection, the build-up of entanglement close to the critical 
value of the field can still be observed in the phase space 
properties. 
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